Decoherence of Schrodinger cat states in a Luttinger liquid 



P. Degiovannj^ and S. PeyssonQ 

Laboratoire de Physique de I'Ecole Normale Superieure de Lyon, 
Unite de mixte de recherche CNRS et Ecole Normale Superieure de Lyon (UMR 5672), 
Groupe de Physique Theorique, 46 allee d'ltalie, 69364 Lyon Cedex 07, France. 

Schrodinger cat states built from quantum superpositions of left or right Luttinger fermions lo- 
I cated at different positions in a spinless Luttinger liquid are considered. Their decoherence rates are 

'n— «' ■ computed within the bosonization approach using as environments the quantum electromagnetic 

' field or two or three dimensional acoustic phonon baths. Emphasis is put on the differences between 

the electromagnetic and acoustic environments. 
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I. INTRODUCTION AND MOTIVATIONS 



Decoherence is a very important issue for mesoscopic systems since it governs the crossover between quantum 
^ ' and quasi-classical transport regimes. Coherence of mesoscopic conductors gives rise to various Aharonov-Bohm 
I ■ interference effects such as permanent currents in a mesoscopic ring, and conductance oscillations as a function of the 
external magnetic field. 

The problem of electron decoherence in metals at zero temperature is an active area of research both from the 
d theoretical and experimentali point of view. Recent experiments claim to observe a saturation ofptke,|dephasing time 
^ T0 at very low temperature*!. Strong discussions among theorisis arose from these observationfflfl'QElEl. The heart of 
I , the debate, summarized fcj£|-esample in Mohanty's recent letterEl, is to determine whether the conventional theory of 
' ^ ' dephasing ia Fermi liquidsBclllJ could explain the saturation of (for .example as an effect of an external microwave 
^ ' radiatiom'EJ) or whereas one should reconsider the theory completelyHjOo. Although such an excitement is a strong 
^ motivation for working on decoherence in disordered mesoscopic conductors, we shall present a model for studying 
I— I the decay of Schrodinger cat states in ID ballistic conductors. 

^-H From a methodological point of view, our line of thought is very close to the one used in atomic physics, for example 

J> in theoretical worksE£lE£l on decoherence experiments in cavity QEDEIL ft relies heavily on the use of simple exactly 
solvable models of decoherence, such as the Caldeira-Leggett modeLJ. We think that this point of view brings a 
different light on the question of decoherence by taking into account the electron fluid as a whole and studying the 
coupling of this many-body system to the external reservoirs. It makes a bridge between atomic physics situations, 
where decoherence is extremely well controlled and for which simplified decoherence models apply almost directly 
' and mesoscopic conductors which are complex interacting systems. In particular, this point of view is well suited to 
] the ID case because of interactions. In a non-Fermi liquid, one cannot keep track of an individual electron because 
^ ■ of orthogonality catastrophe effects. Moreover, electromagnetic or acoustic radiation emitted by one part of the 
Ci ' system and absorbed by another part may have drastic effect on the strongly correlated electron state. These remarks 
^ , motivated our method, based on the study of the coupling of an external environment to the low energy excitations of 
I the ID electron system. Following general ideas of Stern, Aharonov and Imrylij, we have computed the decoherence 
rate of Schrodinger cat states built from localized excitations {e.g Luttinger fermions) introduced at different places in 
the system or at the same place but moving in different directions (left and right moving components). As expected, 
^ ' we have found that such linear superpositions decay into statistical mixtures even at zero temperatures and computed 
• • 1 various decoherence scenarii. This is the main result of this paper and it means that a ID pure ballistic conductor 
exhibits decoherence at absolute zero in the sense of Schrodinger cat states decay. 



One dimensional conductors in the ballistic regime are appropriately described by an effective interacting theory: 
the Luttinger liquicl23. In this effective theory, interactions between electrons are put by means of an electrostatic 
short-ranged potential. Obviously, within the approximation of a linear dispersion relation for particle-hole excitations, 
the Luttinger effective theory contains no source for decoherence. A coupling to an external quantum environment is 
necessary to introduce decoherence in the Luttinger liquid. The quantized electromagnetic field and a 2D or 3D bath of 
longitudinal phonons will be considered in the present paper. The combined use of bosonization and non-equilitucium 
techniques make it possible to solve this problem. Within the bosonization framework, the pioneering workEiJ by 
Martin and Loss investigates the question of equilibrium permanent currents induced by fluctuations of the quantized 
electromagnetic field. They have shown that coupling the Luttinger liquid to QED leads to a renormalization of 
the Luttinger liquid parameters. Our discussion of dynamical and therefore non- equilibrium aspects of the coupled 
Luttinger & QED system will show how this renormalization appears dynamically. Let us mention that coupling a 
Luttinger liquid to one dimensional phonons also renormalizes the Luttinger liquid parameters and drives a ID Fermi 
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liquid to a non Fermi liquid fixed pointc^caCSIia. But a ID phonon bath does not introduce any intrinsic decoherence 
in the Luttinger liquid since, roughly speaking, it does not have enough modes. That's why 2D and 3D phonon baths 
are considered in this paper and we have shown that these baths have enough modes to kill Schrodinger cat states. 

To be more precise on the results presented in this paper, we have shown that the electromagnetic decoherence 
time is much larger, although not infinite, than the natural time associated with the Luttinger system. This is mainly 
due to the weakness and the transversality of the coupling between photons and the Luttinger liquid. Therefore the 
coherent Luttinger liquid paradigm is not ruled out by its coupling to QED. We have also shown that the acoustic 
decoherence is much stronger than the electromagnetic one. This difference comes from the fact that many bosonic 
modes of the Luttinger liquid have an efficient acoustic radiation rate whereas only very few modes dominate the 
decoherence process in the electromagnetic case. In the acoustic case, decoherence takes place over a much shorter 
time scale than dissipation contrarily to the electromagnetic case. 

This paper is organized as follows: the model is presented and the Feynman- Vernon and Keldysh basic tools are 



briefly recalled in section IL Section [II makes contact with the Quantum Brownian Motion problem (QBM). The 



central problem of mutual decoherence of Schrodinger cat states is addressed in section IV. Results are summarized 
and possible extensions are discussed in the conclusion. All along the paper, the electromagnetic and acoustic cases 
are discussed separately. Technical details are gathered in appendices. 



II. ELECTRON SYSTEMS COUPLED TO EXTERNAL RESERVOIRS (QED OR PHONONS) 

A. The Feynman- Vernon-Keldysh method 

Studying an out of equilibrium quantum system boils down to the computation of its density matrix as a function 
of time. A functional integral approach to this problem has been given long time ago in the context of perturbative 
field theory by Keldysht3. Let p{t) denote the density operator for a closed system at time t (in Schrodinger's picture) 
and U[ti,tf] denote the evolution operator between ti and tf, then p{tf) = U[ti,tf]p{ti)U[ti,tf]~^ . The first of these 
operators takes into account what will be called the "forward time branch" and the other one the "backward time 
branch" of the evolution. Evaluating matrix elements of p{tf) can be done with an appropriate propagator which can 
be represented as a double functional integral (one for each evolution operator). Correlation functions are generated 
by introducing external sources coupled to the system's degrees of freedom which have different values on the forward 
and backward time branchea23. Then, one usually takes a trace at time tf. The generating functional obtained this 
way is called the Keldysh generating functional and can be represented as a path integral over a special contour K 
that goes from ti to tf and then back to t^. 

Denoting by ip (respectively ^) fields (respectively external sources) describing the dynamics of the system, and by 
(p-^- and ip- their restrictions to the upper and lower branch of Keldysh's contour (respectively ^_), we have: 

where p{ti)[ip^\t- , </3-|tJ denotes the kernel of the initial density operator. This is the basis for Keldysh's non equilib- 
rium perturbation theory. Because of the doubling of degrees of freedom, one obtains a 2 by 2 matrix Green function 
which can be related to well known Green functions by: 

Gi.,y) = (a.(^,.)),,.).i.-r> = ( cfilll] ) 

Here T denotes, the usual time ordering, T the anti-chronological time ordering, < and > the lesser and upper 
time orderingfS. For bosonic oscillators initially in thermal equilibrium, Green's functions are explicitly known (see 
appendix ^ . 

For a matter system coupled to QED, we are interested in the evolution of the reduced density matrix for the 
matter system. Bosonization techniques enable us to treat the electron fluid through an effective bosonic theory. Our 
strategy will then be to integrate over the electromagnetic fleld degrees of freedom and to deal with matter degrees of 
freedom in a second time. It is convenient to choose the Coulomb gauge, in which dynamical degrees of freedom of the 
quantum electromagnetic field (transverse photons) are decoupled from the Coulomb interaction. The latter is taken 
into account through the effective action for the matter system which uses only instantaneous interactions such as the 
Luttinger liquid description. Integrating out transverse, photons gives a non-local functional integral of the current 
density, called a Feynman- Vernon influence functional. For the sake of simplicity, we shall assume that these two 
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systems are independent at initial time (the initial density operator factorizes p{ti) = pMatter Pqed) and that the 
electromagnetic field is initially at equilibrium at temperature T = l/ksP- In the present case, using bold symbols 
for vectors and matrices, the influence functional is nothing but the influence functional for the current density: 



where is the appropriate functional kernel representing the density operator for the quantum electromagnetic 
field at given temperature. The functional integral ^ can be evaluated in terms of Keldysh's Green function for the 
electromagnetic field at finite temperature D'^f!{x,y) = —i{TKA"{x)A^,{y)) and is equal to: 

^b+,J-]=exp|-^ ^ ee' fd^xd''yj:{x)D:f!{x,y)/Ay)] (3) 

V (e.e')e{l,-l} / 

Noise and dissipation kernels are defined by: 



v"P(x,y) = --^{Dfix,y)) (4) 
1 



V"'i^,y) = --{D-^\x,y)-D'X^{x,y)) (5) 



^[j+J-] = exp J ( 



and we have: 

d^xd^y'{j+ix)-j_ix)).,^{x,y).ij+iy)-j_iy))^ (6) 

X exp{^iJd'^xd'^y\jAx)-j_{x)).r]{x,y).{jAy)~j_{y))^ (7) 
Correlation functions can be obtained from a generating functional which takes the environment into account: 

Within the elastic approximation, coupling to acoustic phonons can be treated along the same line. Longitudinal 
phonons create a potential V{x,t) such that eV{x,t) = Dp{x,t)div{u{x,t)) where u{x,t) is the elastic deformation 
field, p{x, t) is the electric charge density and D denotes the electron-phonon coupling energy. Longitudinal phonon 
dynamics is described by a quadratic action: 

S[u{x, t)] = ^ J {{dtu)^ - c| (div(M))2) d^xdt (8) 

where cs denotes the sound velocity and pM the volumic mass. The resulting phonon influence functional is a Gaussian 
in terms of p±{x) with kernel D'^'^ {x,y) = —i{u"{x)u^,{y))i3: 



^[p+,P-]^exp\~^-^ f d'^xd'^y ^ Peix)pe'iy)id^,ady^f3V"-'^^){x,y) \ (9) 

(e,e')6{ + ,-P 

In a generic electron sjjstem, this influence functional a priori contains quartic fermion terms. The one, dimensional 
interacting electron gasE^I is quite interesting since it can be described by an effective free bosonic theor]fiB. Neglecting 
environment-induced umklapp processes, the influence functional is Gaussian (although non local) in terms of the 
bosonic field. Such terms arise from the introduction of 2/cf components of the charge or current densities in the 
environment's infiuence functional. There is a factor exp {2ikp{cr + a')) in front of the term involving two tp]^tpLS. At 
incommensurate filling, the usual averaging argument can be invoked to rule out these Umklapp contributions. Other 
terms contain both tpli'^'L and 4''l4'R and would imply momentum transfer of order 2kp (except for special values of 
filling). In QED's case, 2kp photons would have a very high energy compared to the Luttinger typical energy hvs/L. 
In the acoustic case, since kpL >> vs/cs, it would also be the case. That's why we will not take into account 2kF 
components of charge and current densities. Thanks to a Fourier mode analysis, this problem boils down to a set of 
independent harmonic oscillators, each of them being linearly coupled to a bath of quantum osciUators. This problem 
is known under the name of Quantum Brownian Motion (QBM) and has widely been studiedEjEil. Necessary results 



will be recalled in section III. Elementary excitations of the Luttinger liquid can be created using vertex operators. 
States created this way are nothing but coherent states. Since the evolution of coherent states in QBM can be exactly 
computed^, so can the evolution of elementary excitations in the Luttinger liquid. 
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B. Mode decomposition for the Luttinger liquid 

Let us consider a Luttinger liquid on a circle of length L, with both left and right chiralitics. As inEl, a denotes 
the coordinate along the circle. Low energy excitations of the Luttinger liquid are described by the theory of a free 
compactified boson, whose action is given by 

S[^] ^^JdcTdt (vs'idtipf - vsid,^)') (10) 

The field is compactified on a circle of radius Rc, and the Euclidean path integral representing the finite temperature 
equilibrium partition function contains a topological term taking into account zero mode quantization. The space of 

states of this free bosonic theory is a representation of a U{1)^ x algebra (two commuting copies of an oscillator 

algebra): [Jk,Ji] = fc<5fc,-;l and [Jk,Ji] = kSk^-il- In Wen's workS, J; and J; for I ^ are called hydrodynamic 
modes, whereas the Jq and Jq are called zero modes. In term of these modes, the Hamiltonian of the Luttinger system 
is given by: 

Htot = ^Y.^JiJ-i + JiJ^i). (11) 

lei 

The finite temperature equilibrium partition function with a magnetic flux x^/^ is given by: 

J_ , „_ -Ib", 

l^(r)| 



Lutt(x) = ^TZT^ 2^ (12) 



a)e(Z/2)xZ 
m (mod 2) 



where a = gR^ encodes the interaction strength {a — 1 for free electrons) and 

Pn.m^ny^+—= p,^„^^n^/^-—=. (13) 

States of interest will be U{1) x U{1) highest weight states: with Jq (resp. Jo) eigenvalues Pn,m (resp. Pn,m)- The 
charge and current densities can be expressed in terms of the U{1) modes as follows: 

p{a,t) = ^(j„(i)e2-r- + j„(i)e-2-f-) (14) 



ne'i 

j{a,t) = £^^(j„(t)e2-f--J„(t)e-2..j.) 



(15) 



The jn, to)^^, — \pn,m+2xT Pn,m+2x) State carries a charge 2n G Z and a current (m + 2%) us/LoHH. If we assume that 
our Luttinger system is realized using the edge states of a FQH fluid, then the two edges of the sample are polarized 
with charges {a = l/iy is an odd integer): 

(771 \ / 777 \ 

77 + — j and = e (^77 - — j . (16) 

Obviously, U{1) descendants carry the same global charges and current than the highest weight states. Remember 

that JiXU-iprimaries are created using vertex operators Vn,m{o',t) which are normal ordered exponentials of bosonic 
modea^j'tZl. To complete this brief description, let us recall that the original fermionic operators get renormalized 
(orthogonality catastrophe historically introduced by AndersonEa) , and that the corresponding renormalized fields are 
the so-called Luttinger "fermions" which correspond to 77 = ±1/2 and 777 = ±1. Within the context of fractional 
quantum Hall effect on a cylinder - or an annulus y-e-dge fermions carrying a charge localized on one of the two edges, 
can be created or destroyed using vertex operatora^J with 77 = ±1/2 and 777 = ±v~^ . Typical experiments which may 
be performed on edge excitations of an annular 2DEG in an AsGaAl heterostructure lead to the following numerical 
values which we shall use in the rest of this paper : vs —i- lO^ms^^, L ~ 30/im, a = 3 and T ~ — 10 iiT. 



C. Coupling of modes to tlie environments 



In order to respect cylindrical geometry, the environment field is supposed to live in a cylindrical cavity with the 
same revolution axis than the Luttinger circle. The cavity is of radius R and of height h and V = ■nR-'h will denote 
its volume. Periodic boundary conditions are assumed in this vertical direction. 
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Electromagnetic modes are quantized using their momentum along z axis kz^ their angular momentum moment 
around z axis I G Z and the number n of radial zeroes of the electric field. The boundary condition on the cavity 
plays the role of polarization and accounts for the helicity degeneracy. Transverse electric (TE') modes have their 
(transverse) electric field orthogonal to the cavity's edge, and transverse magnetic {TM) modes have their magnetic 
field orthogonal to the cavity's edge. Normalized expressions for the modes can be found in any book on optical fibers 
and are recalled in appendix |^. Of course, the box is considered as large meaning that electromagnetic modes form 
a continuum compared to the Luttinger modes. 

Phonons are quantized according to the vanishing of displacement field u on the cavity's boundary. Longitudinal 
modes can be written as a gradient u = ^(ip). As shown in appendix longitudinal modes are classified by (fcz, I) 
and n where n is the number of zeroes of the ip potential. 

It is convenient to introduce some kind of "phase space coordinates" to describe the I > 1 modes of the Luttinger 
Liquid : 

ql^^ ^ ^^{Ji + J-i - Ji - J-i) and = -i=( J; - J_; - J, + J_;) (17) 

-iql'^ = ^^{Ji ~ J-i + Ji ~ J-i) and = -L(J, + J_, + J; + J_;) (18) 

Formulating the problem in this way shows that the interaction Hamiltonian is nothing but a linear coupling between 
harmonic oscillators (qx coupling). More precisely, for each Z > 1, two harmonic oscillators are coupled to the acoustic 
or electromagnetic modes of angular momentum I. The two Luttinger modes of fixed I are coupled to different sets of 
environment modes. The same conclusion is true for zero modes: the global current couples to ^ = electromagnetic 
modes whereas the global charge decouples from the field's propagating modes (Coulomb gauge effect). In the acoustic 
case, only the total electric charge couples to Z = acoustic modes. 

The influence functional for each of these environments can explicitly be obtained in a nice form. First of all, the 
charge and current densities can be expressed as linear combinations of the above phase space coordinates. Then, 
performing integration over the vector potential and spatial coordinates leads to the QED functional (assuming that 

= qo and g^' - 0): 

^PV = cxp(-|^5] Y: ee'^Vjk f dt ds qlf^t) G'rj\t - ^) llf' i-'^) \ (19) 



i , >0 



(e,e',Q)e{ + ,-}3 



where the dimensionless coupling constant only depends on the fine structure constant uqed, the Luttinger liquid 
interaction parameter a and the ratio of velocities vs/c: 

g = ^^^QER,r_lY (20) 

a \ c J 

This mode expansion is related to a multipolar expansion. Luttinger modes of momentum I contribute to the electric 
(resp. magnetic) multipolar expansion starting at order 21 (resp. 21 + 1) In particular, / = is mainly a magnetic 
dipole and an electric monopole (corresponding to the global charge and current around the Luttinger ring). The 
I = 1 modes contribute first to the electric dipole and the magnetic quadrupole. 

Longitudinal phonons can be treated in the same way. In this case, one needs to perform a canonical transformation 
on the Luttinger phase space coordinates, exchanging and Q;'^' in order to obtain a Gaussian expression similar 
to |l^. The dimensionless coupling constant corresponding to longitudinal phonons in dimension d = 2, 3 is given by: 

where cs is the sound velocity, pM the volumic {d — 3) or surfacic {d = 2) mass and D the typical electron/phonon 
coupling energy. With typical values D ~ 7 eV , cs — 3. 10'^ ms~^ and pM — 3. 10'^ kgm~^ , we get, in dimension 
three gph{L) ~ 4. 10"^. 



6 



III. RELATION TO THE QUANTUM BROWNIAN MOTION 
A. The Quantum Brownian Motion model 

This basic model consists in a single quantum harmonic oscillator (called "the system") coupled to a bath of 
oscillators (called the "reservoir" ) : 

L^^ie- nV) +J:Ciiii + T.^ ^iii) (22) 

We assume that these two systems are initially independent and that the oscillator bath is at equilibrium with inverse 
temperature /3. One could also assume that the whole system is initially at equilibrium but it rftakes things more 
involved without illmninating the discussiorL—. Interested readers are refered to Grabert's reviewKI, to Hakim and 
Ambegaokar's paperEJ or to Romero and Pazo. The question addressed here is to compute the time evolution of the 
system. The physics depends on the reservoir's influence at each frequency, encoded in appropriate spectral densities. 
Therefore, several natural time scales will appear in the Luttinger & environment problem: 

• The cut-off time : above a certain UV cutoff, the effective description by a Luttinger liquid is no longer valid. 
Another natural cutoff could also be provided by the spectral distribution itself (Debye frequency for phonons) . 
For this reason, an UV cutoff is present in the model. The associated time will be denoted by 1/A. Below this 
time scale, our simple model cannot be considered as valid. 

• The environment time scale te is the time needed by the light {tem — L/c) or by phonons {ts — L/cs) to 
circle around the system. 

Spectral densities and influence functionals can be normalized with respect to this characteristic time (they 
only depend on ujte)- The low frequency regime is defined by lote << 1 and the high frequency regime by 
ute » 1- 

• The Luttinger time scale tl = L/vs is the time needed by one excitation of the Luttinger liquid to circle 
around the system. This is the natural time scale from the Luttinger liquid point of view. 

• The thermal time scale tth = h/ksT is the inverse frequency associated with temperature T. The cor- 
responding thermal length is l{j3) = CTTh for photons and = cstth for phonons. Let us notice that at 
T ~ IK, l{p) = 1 mm for photons whereas it is of order 10 fiia for phonons. 

The general solution to the Quantum Brownian Motion problem is due to Hu, Paz and Zhang_(HPZ)ii who 
computed the evolution kernel for the system's reduced density matrix. Caldeira and Leggett's workEfl is concerned 
with a special case of reservoir, corresponding to the so-called Ohmic spectral density J'{uj) oc cu. As we have already 
noticed, the Luttinger liquid can be seen as a collection of harmonic oscillators and some zero modes. As we shall 
see, each harmonic mode coupled to the environment is a QBM problem. Each reservoir is characterized by a set of 
spectral functions (one for each I >0), which are made dimensionless for simplicity. 



In the electromagnetic case, using expressions p36| and B7 in appendix M, we obtain 



Jii^)^9^ E 0^'^' 6{tem{uj-u;i)) (23) 

y , 77^ , 2. U! ITEM 
I I ll=l 



The I — case needs a slight modification: 



L3 ^ 2)| 



M^) =9tT 2^ KrEMico - LOi)) (24) 



The electromagnetic influence functional can be rewritten as: 



*IlE E / dtJ^dsiql'^^+-ql'^^-mr,,{t-s)iql-^++ql"^-)is) 

/ 2 oo rtf pt 

-liE E / '^^Z dsi^l 



xe^p -72E E dt dsiqr^-qr-mMt--s){qr-qrnis) (25) 
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Dimensionless dissipation and noise kernels for each mode are expressed in terms of the spectral density by: 

rji{s) = / (Lj Ji{uj)sui{llis) (26) 



viis) — — duj Jiiuj) coth ( — — I cosfws) (27) 
c Jo V 2 / 



An acoustic reservoir leads to similar expressions for spectral functions (using expressions CI and C2) 



Ji{cj) = gp,^{L) ^ J2 liM]{Tsu:if5{Ts{u:^uoi)) (28) 

// l!=l 



The I = case needs a slight modification: 



M^)=gph{L)^ MUrsu:if5{Ts{u:-uJi)) (29) 



Numerical computations of spectral densities as well as analytic estimates of their asymptotics are available (see 
appendix The electromagnetic case is illustrated on figure 3. In this case, alH 1 modes are supraohmic {Jiioj) 
decreases faster than w) at low frequency. The / = 1 modes sh ow a n ohmic behavior (J7i(a;) oc w). In this case, the 



dissipation kernel is local in time and, as we shall see in section IV F , an effective Caldeira-Leggett model can be used 
to perform analytic computations. 

In the acoustic case Ji{Ld) goes as [lotsY'^^'^ in the low frequency regime and as {u)tsY~^ at higher frequencies. The 
main difference between the electromagnetic and acoustic reservoirs is that the natural Luttinger frequency 2'it/tl 
falls in the low frequency domain for QED, whereas it does not for phonons since ts/tl >> 1. 

B. Phase space evolution using Wigner functions 

The time evolution can be computed using the Wigner function associated with the system's density operator. This 
form is especially adapted to the study of decoherence of Gaussian wave packets. Moreover, it provides a nice quasi- 
classical insight on the evolution of the system since in the Luttinger liquid, charge and current density fiuctuations 
play the role of "phase space coordinates" for the hydrodynamic modes and within the FQH effect framework, encode 
the shape of the incompressible quantum Hall fluid dropleta. The Wigner function associated with an operator B is 
defined by: 

WB{p,q)= J dye'Py/^{q~^\B\q+^) (30) 

We use the following notation for phase space: (f) = ^ ^ ^ • evolution kernel for the Wigner function can be 

computed from the density operator evolution kernel and is given by: 

Jw{p, q, t\po, go, 0) = M{t) exp J*(^- + ^o- </>o)- A{t). {N. cj, + NqM^ (31) 

where the Nq and N time dependent matrices are given by: 

The Ui{t) {i G {1,2}) functions are defined in Hu, Paz and Zhang's paperEl as solutions to the classical equations 
of motion with dissipation and boundary conditions (ui(0), ui(i)) — (1,0) and (m2(0), M2(i)) — (0,1). The A{t) — 
(ai_j(t))(ij)g{i_2}2 matrix is defined by 



.]{t) = \ [ dsi j ds2Ui{si) visi - S2)Uj{s2). 

^ Jo Jo 
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The reader can check that, when the reservoir decouples, the evolution kernel reduces to a delta function, giving back 
the classical evolution in phase space. Henceforth, we shall call —N^^Nq the "Wigner evolution operator" and denote 
it by Ut- When turning on the coupling to the environment, the classical delta distribution in phase space spreads 
into a Gaussian, the center of which moves according to Ut- As we shall now see, the evolution kernel encodes all 
effects of dissipation and decoherence. 

The evolution of a density operator built from Gaussian wave packets can be computed exactly. The generic form 
of a Gaussian Wigner function is, up to some normalization: 

W[4>, K, Q] (</.) = exp (^-i*(<^ - 0). Q. (0 - ^) + z*K. </.^ (33) 

The 2 by 2 matrix Q encodes the spreading of the packet, (j) is the center of the packet and K a phase modulation. 
Gaussian wave packets and in particular coherent states lead to Gaussian Wigner functions. A Gaussian Wigner 
function remains Gaussian. If W{t = 0) = W[(f>o, Ko, Qo] then, after time t, 

W{t) = e-'^[^°'«°l(*) X W[Ut. <po, K{t), Q{t)] (34) 

where the basic parameters at time t are given by d[Ko, Qo\{t) = *Ko. D{t). Kq and: 

D{t) = {2Qo + 'NoA{t)-^No)~^ (35) 
K{t) = %-\{l-2QoD{t)).Ko (36) 
Q{t) = %-KiQo-2QoD{t)Qo).Ui' (37) 

The D{t) matrix contains the decoherence effect. To understand this, let us start with a coherent superposition of 
two Gaussian wave packets {tpi + ip2)/y/2. The initial Wigner function is given by: 



W{0) = i (1^11 + Wi2 + W21 + W22) where W„^(0) = W[4>o,, ^p, Qo] = W 



2 

Under time evolution, the form of the wave packet is preserved up to a global factor and a phase modulation: 

Wa^it) = e-<^-^(*)+^®"^(*.« W[Ut. 4>a,Ut. 4>0, Q{t)] (38) 



where: 



dap{t) = d[-ia\{(l>a-(l>p),Qo\ (39) 
Qa0{tA) = 'Kgit).{4>-Mt)) (40) 



K,it) = Kit)-p^%~\Ko (41) 



Each density operator \'ipa)(jpa\ has its own evolution, described by (Pait) = Ui.(j)a and Q{t). Note that dii(t) = 
d22{t) = 0. The coherence part is contained in W12, and evolves according to Ut.(j)i2 and Q{t) plus an exponential 
factor e-'^i2(t)+iei2(t,<^)_ T^e di2{t) factor gives the attenuation of off-diagonal correlations and should therefore be 
interpreted as the decoherence factor between the two wave packets. 

Unfortunately, as noticed before, explicit and closed expressions for Ut and D{t) are not known for a general supra- 
ohmic environment. Therefore, we shall perform a perturbative expansion in the coupling constant. In the case of 
the decoherence factor di2{t), it is enough to start from standard harmonic oscillators expressions for the Ui and U2 
functions since A{t) is of first order in the coupling constant. The decoherence matrix can then be expressed as an 
integral over the spectral density: 

D{t) = r J{uj)D{n,uj,t) coth duj. (42) 



where D{Q.,ijj,t) has an explicit but involved expression. Let us discuss its asymptotics in various physically relevant 
limits: 
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• At short times fit << 1, can be neglected provided the spectral density of the bath contains modes with 
frequencies much higher than f2. Then, this regime is dominated by the high frequencies of the bath and: 



D{n,u;,t)^{ 2 i (43) 




The last expression is only valid at very short times and contains non-Mar kovian effects. Even in the Caldeira- 
Leggett model, at very short times, all frequencies of the bath take part in the evolution of the system, leading 
to these non-Markovian effects. 

• When the condition Qt << 1 is no longer valid, dissipation effects with exponential relaxation will appear 
through linear terms in t. We must of course assume that t is much smaller than typical relaxation times. At 
growing times, linear terms will dominate oscillating resonant ones and provide a "Golden rule" estimate for 
the relaxation and decoherence times. This approximation will be used in the next section in order to evaluate 
the typical decoherence time of Schrodinger cat states in the Luttinger model. 

When t reaches the dissipation time scale, this perturbative treatment breaks down. The long time regime of 
decoherence can however be computed in some cases using the Caldeira-Leggett model. 



IV. MUTUAL DECOHERENCE OF ELEMENTARY EXCITATIONS IN A LUTTINGER LIQUID 

A. Statement of the problem 



As recalled in section IIB, elementary excitations of the Luttinger liquid are created by vertex operators Vn^mi'^)- 



Such a state is not an eigenstate of the Hamiltonian but rather corresponds to the introduction of a "localized" 
excitation at point a around the circle. Let us consider a Schrodinger cat state build as a superposition of a right 
moving Luttinger fermion (n = 1/2 and to 1) at different places around the circle: 



In an isolated system, such a state will evolve according to: 

1 



l^flfl(t)) = {i^l,iai,t)\0)+iPl,{a2,t)\0) 



(44) 



(45) 



Therefore, this coherent superposition will remain coherent, as it should in any isolated quantum system. Switching 
on the coup liri j g t Oj thp quantum electromagnetic field or phonons changes the situation: according to general works on 
decoherencccJcJcj'cS, we expect this Schrodinger cat to decohere into a statistical mixing of two states: one excitation 
at one position, or the excitation at the other position. Here, two questions will be addressed: what is the strength 
of the decoherence process and on which time scale does it take place? 

In the following, two cases will be considered: the R/R Schrodinger cat, already presented in equation |44| and the 
R/L Schrodinger cat, defined as: 



|V«L(0)) = -i=(vlj(ai)|0)+^I(a2)|0) 



(46) 



Practically, the case of zero modes is simpler and will be considered in the next paragraph. We shall then turn to the 
I 7^ modes in se ction [VC. Decoherence time estimates for electromagnetic {resp. acoustic) reservoirs are given in 
IV D {r esp. VE ). An effective Caldeira-Leggett model will be used to deal with the long time behavior in sections 
iVF toffVH. 



B. Evolution of zero modes 



Eor the zero modes, the zeifo, coupling evolution corresponds to a free particle (and not to an harmonic oscillator). 
The Ambegaokar and Hakimcll method (exact diagonalization for the coupled system) can be used to compute the 
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evolution of the density matrix for the zero modes. The corresponding exphcit formula in the Luttinger liquid case 
is, at finite temperature: 

(n, m\p{t)\n', m!) = {n, m\p{0)\n', m') x exp {-i{ujn,m{t) - ujn,rn{t)) t) x exp (^--^(^" - ™')^^ (47) 

where hjJn,m{Q) denotes the energy of \n,m) in the isolated Luttinger system and: 

TTwc o /•+00 f, , , / sinfo^i) \ duj , , 



2La Jq nvs \ ivt ) ijO 

f^°^ . ^, . , //3fiw\l-cosM 
d{t) = / dujMu) coth ^ I > 49 

TEM Jo V 2 / (Ujt)^ 

The first term in ^ can be interpreted as a dynamical renormalization of vs/a. Only vs/a is renormalized since the 
charge density does not couple to the transverse degrees of freedom of the electromagnetic field. To be precise, at 
t — > +O0, the velocity and interacting parameters of the Luttinger liquid get renormalized as v'g = vs- C and a' — a/<^ 
where the renormalization constant C is equal to the t +00 limit of: 

at)' = 1 + — / jo(^) - 1 — (50) 

■KVS Jo \ Ujt J LO 

The renormalization effect is of course strongly cut-off dependent. The dimensionless coupling constant appearing 
here is ^ ~ IQ-^. 

Q C 

The second term is the decoherence coefficient between two different highest weight states |n, m) and \n' ,m') of 
the LL. Decoherence takes place in a time of the order of the cutoff time and then reaches saturation. Figure 4 
summarizes d{t) and C(^)'s behavior. 

The typical value of the t — > +00 value of the decoherence exponent is of typical order g: 



+00 



f I3huj 
Lw2 """" \^ 



I _ Jo(a;) coth ( ^ ) (51) 



In the acoustic case, computations can be performed in the same way. Since acoustic zero modes couple to 
the total charge of the Luttinger system, the decoherence factor between states \n,m) and \Ti',m') is found to be 
exp {—2ad{t){n' — in'Y)) where d{t) is obtained from ^ by using the acoustic spectral density and rg insted of their 
electromagnetic counterparts. The Luttinger parameters a and vs also get renormalized. In the acoustic case, only 
vsa is renormalized. Then, a' — aC,ph and v'g — vsCph where Cvh is the t +00 limit of Cphit), obtained by using the 



acoustic spectral density and the speed of sound in formula 50 



In both cases, the final decoherence exponent is proportional to the square of the difference between the total 
current (resp. charge), quantities which measure the "distance" between the two quantum states. Such a result is 
expected since, as explained in C. Cohen Tannoudji's lecturescll, such a dependence is common in the case of a linear 
coupling with a conserved quantity. In particular, a Schrodinger cat obtained by superposing the same elementary 
excitation of the Luttinger liquid at two different positions along the ring has all its decoherence due to hydrodynamic 
modes ! We also notive that zero mode decoherence has a weak dependence in the cut off and temperature. 



C. Spatial dependence of decoherence 

Using the explicit form of vertex operators, one easily finds the relevant parameters to be used for the decoherence 
of each mode. Of course, these parameters depend on positions of each of vertex operator. Here, we shall only present 
the results for R/R and R/L Schrodinger cats (cri2 = (Ti — (72)- 



+00 

(52) 



dnLit,a„a2) = igi {m^<(t) (l - sin^ (^) 



(53) 



+ An'a'Di'^it) sin' + 2nmaD['Ut) sin f^^) ^ (54) 
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Here D^^\t) denotes the decoherence matrix for the q^^^ modes computed along the hnes of section [II B. The main 



change from the HPZ computations arises from our normahzation choice for the g/S. The effective spectral density to 
be used in the HPZ formulas is given by: 

c 

This rescahng takes into account the ratio of the Luttinger mode I time scale {i.e. ti^/1) and of the environment time 
scale TE = L/c. 

The appearance of an odd dependence - in term of ai2 - in the d_RL(i, ci, (72) coefficient is understood by noticing 
that an appropriate parity operation transforms the R/L Schrodinger cat into an L/R one. Therefore cti, (T2) 

is invariant into simultaneous changes cti <-> (T2 and nm i-^ —nm. 

A first estimate is obtained using a perturbative approach. We perform a secular approximation and retain only 
terms linear in time. The corresponding decoherence rates are given by: 



T = Stt An.m ■ sm — - — 



(55) 



T/^/^^(ai,a2) V m2+4n2a2 \ L JJ ujite 

Here A„^m is the conformal dimension of the vertex operator Vn^m{o'). Not surprisingly, the decoherence time of a 
R/R Schrodinger cat diverges when 0-12 0. This result is obvious since in this limit, the initial state is a pure state. 
For L/R cats, the decoherence time shows a slow variation in term of the differences of positions. 



D. Decoherence time estimations: QED's case 

Using asymptotics of spectral densities (see appendix ^ , one obtains the decoherence time of the Ith modes in the 
R/R case: 



f2TTVsY'-^~^^ + .2/7r?0-l2\ 



Similarly, for the R/L case: 



= 45A„,™ ^ 1+ ^2 , .„2„2 CoM^^ (58) 



T/'*/'^)(ai,a2) " "'""K c J ■(2/ + l)!-V + Aa^n 

Since vs/c ~ 10"'^, decoherence times for the I and / + 1 modes are related by a typical factor of 10^. This argument 
shows that the / = 1 modes dominate the decoherence process. Physically, higher Luttinger modes contribute to 
higher electric and magnetic multipoles, for which radiative dissipation is known to be weaker. Since dissipation 
governs decoherence, this is the physical reason for the predominance of the I — 1 Luttinger mode in the decoherence 
process. The decoherence time of the I = 1 mode is nothing but the electromagnetic relaxation time: 



16- OQED /«5y ^ ^0-8 1 

3 a V c / 



Numerical results for the decoherence times are shown on figure 5 for Luttinger fermions. 

The temperature dependence can be found easily since coth ((3hLu/2) varies slowly around uji in a scale gcui. There- 
fore: 

^'"'"'^'^^'""'^^ -tanhf^^) (59) 



r/''/''^(ai,a2,T = 0) V2fcsT 
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E. Decoherence time estimations: acoustic case 



Since the sound velocity cs is much smaher than vs, the condition loite << 1 docs not hold for the coupling to 
phonons. Explicit computations show that the behavior of spectral densities for the coupling to longitudinal phonons 
differs from the electronic one. Luttinger modes frequencies fall into a range of frequencies where the acoustic spectral 
densities are proportional to (u!tsY~^ ■ Remember alsa,that the natural cut-off frequency for the phonon bath is given 
hy lod — Cs / a where a is a typical microscopic lengtfO. The zero temperature acoustic decoherence rate of Luttinger 

modes is given by formulas |5^ where the r|'*'' = J'i{uji)/ujits factor is given by: 

p(3) ^ 9ph{L)P vs j-gQ^ 
' 2 Cs 

rf) = (61) 

In opposition to QED's case, these damping rates do not decrease with increasing I. In the QED case, higher I modes 
are bad antennas for the microwave radiation emitted by the system. In the case of phonons the situation goes the 
other way because of the longitudinal coupling. Although in some cases the coupling constant (7p^(L) is very small, 
"decoherence repartition" effects between modes plays a much more important role here than in QED's case since 
one has to sum up over many mode contributions to decoherence. For Luttinger modes of index I, the perturbative 
expansion is governed by the relative damping rate 7; /w; , an upper value of which is given by 

M^^^^ ford = 3 and ^I^Jord = 2 (62) 

where Qphia) is the rescaled coupling constant for the length a (typically of order 1). Assuming that ujdts — L/a is 
much greater than one, we see that all I ^ Q modes can he considered as weakly damped. 

The total decoherence exponent in the linear regime is obtained by summing over all the modes up to the Debye 
frequency. For R/R Schrodinger cat states, one finds: 

r,.^^^'^\a,,a,X) = 8,r^A„,„ g pj'^) coth sin^ (63) 

Since we sum over a large number of modes, the decoherence time rapidly decreases when (J12 » avs /cs, a spectacular 
effect due to the ratio cs/vs << 1. Roughly speaking, the Luttinger fermion has the time to circle many times around 
the loop before emitted phonons escape whereas it barely has the time to move in the electromagnetic case. This 
"averaging effect" explains why the dependence in the initial relative position is much weaker for acoustic than for 
electromagnetic decoherence. The maximal inverse decoherence rate can be expressed as an integral in the limit 
L » avs/cs {kB&D = hujoY 

The temperature dependence is very weak (remember we are typically working in situations where lodtl > 10^ and 
UbT ~ hiUL)- To be precise, it goes like: 



^Wfi)(r)-7W^)(o) ^ /j^v" 

In opposition with the photon bath case, the total acoustic decoherence time scales as L^^ in units of tl 



(65) 



F. Caldeira-Leggett computations 

The Caldeira-Leggett model corresponds to the Ohmic spectral density: at low frequencies, J{lo) = ^'^J" . In this 
case, for time scales large compared to the cut-off time, noise and dissipation kernels are local in time. The solution 
of the model is then much simpler. In the HPZ approach, the equation of motion defining the {ui)i^i^2 functions 
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can be solved exactly, taking into account non perturbatively all effects of dissipation. For general spectral densities, 
solutions of the equation of motions are given in full generality by a Laplace transform of the form: 



i[p) = ^ , and rifp) = / dto^ ^. 

a 

inverse Laplace transformer: 



This expression has clearly a cut along the p-imaginary axis for \p\ < ojc (^c is an UV cutoff such that J^{lo) = 
for uj > LUc). It has no poles in the physical sheet. The Bromwich contour that encircles the cut is used to find the 



Siw] 

The self-energy due to the bath is given by 



u{t) — 2 j S{u!) sm{ujt) duj 







^R{Lu) + i^i{Lu)^2FP f "^/^"^i diu' + in sgniu) J {\lo\). 



(66) 



In the weak coupling regime {"fTL << 1), a Breit-Wigner approximation can be performed since ^/(cj/j) << w^. 
Within this approximation, the m functions can be approximated by standard damped oscillator solutions. In fact, 
one could then imagine to use an effective Caldeira-Leggett model in order to estimate the decoherence properties of 
all Luttinger modes. 

In the electromagnetic case, the decoherence process for the I > 1 modes takes place on a typical time scale of order 
T; ~ Ti(c/ws)^'~^ which is much longer than for the Ohmic I = 1 modes. In our particular problem, using asymptotics 



D7, electromagnetic dissipation for the dominant I — 1 modes is given by 7 = gr^^ /3. An effective Caldeira-Leggett 
model for this mode can then be used. In the acoustic case, for a two dimensional phonon bath and within the relevant 
range of frequencies lo^ < uj < lojj, the spectral density J^i{uj) is ohmic for any I, the Caldeira-Leggett model can be 
used to describe the long time dependence of decoherence. 

Although the HPZ method can be used to perform explicit computations, an approximate master equation can 
often be used to obtain the time evolution of the density matrix of a small subsystem S coupled to a "reservoir" TZ. 
The usual approximation consists first in assuming that the state of the reservoir is unaltered by the coupling to the 
system, and to forget correlations between S at time t' < t and TZ at time t. The second approximation usually done 



consists in neglecting non-Markovian terms which can occur for example in the noise keri, 
always local in the CL model). The condition for these two approximations to be valid isl 



(the dissipation kernel is 



T^. Tr (hIt^ psnm (67) 

where Hsn denotes the coupling between the system and the reservoir, psn{(i) the initial total density operator and 
Tc the correlation time of the reservoir. In the present context, this condition can be expressed as: 

6 Jiix/TE)coth (^-^x) dx « (^] (68) 



2L J \tcJ 

where is a dimensionless coefficient that characterizes the spreading of the l-th mode initial state : = 1 for any 
coherent state and = coih {jrhlv $ P / L) for thermal equilibrium. 

Within the temperature range used here T/Tl ~ — 10 {Tl — 2Tihvs /kBL), we have /(/3) >> L for the electro- 
magnetic bath. Therefore, an estimate of the correlation time for the QEM field is given by the thermal time. High 



frequency asymptotics D£ show that temperature dependence has negligible effect on the l.h.s. of |6^ (the effect could 
be important only for the / — 1 mode). Then, within our temperature range, the l.h.s is approximately bounded 
by gL/a. The r.h.s is bounded by L/a at high temperatures and (vs/c) {T/Tl) at lower temperatures. Therefore, 
since g ~ 10^^, the validity condition is valid for the photon bath at temperatures above g {L/a) (c/vs) x T^. In the 
acoustic case, the l.h.s. temperature dcpendance is only relevant for values of a; below (cs/vs) x [T /Ti). But the main 
contribution to the l.h.s comes from higher frequencies is of order gph{L) {L/aY — gph{<i) {L/a). At low temperatures, 
the r.h.s can be bounded by {cs/vs) {T/Tl) using h(3cs as an upper value for Tc- Condition ^ can therefore be 
rewritten as {Lvs / csa) gphio-) << T/T^. Subsequent computations will assume that the acoustic UV cutoff is much 
higher than the typical Luttinger liquid frequency 2'kvs/ L, i.e that Lvg/csa >> 1. The master equation can therefore 
only be used for quite high temperatures. 
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In the case of an underdamped oscillator, as pointed out by C. Cohen- TannoudjiO, the approximate Markovian 
description can be used provided the noise kernel is local compared to the observation time scale. In the present 
case, this means 7 << luq and ksT >> hj. In the limit of very low temperatures ksT << Hj, this description is 
still valid but one could expect very long time algebraic tails which precisely take into account non Markovian effects 
induced by the divergence of the bath's correlation time ft/3. This possibility is discussed in appendix ^ where these 
corrections are computed in a weakly damped regime and shown to be extremely weak. 



G. Long time decoherence (Markovian master equation treatment) 



Within the secular approximation, the master equation can be translated in the following evolution equation for 
the characteristic function. In the case of a single oscillator of rcnormalized eigenfrequency f2 coupled to an Ohmic 
bath with dissipation rate 7: Zt{X, A) = TT{ps{t) e'^°* e^°) is given by: 



2) dx 



Zt{X,X) 



7 



ZtiX.X) 



The general solution of this equation is given by: 



Zt{X, X) = Zo{ Ae-*+*"S Ae^*-'"*) . exp 



1 - e"^* 



XX 



(69) 



(70) 



Such a formula immediately shows that, at zero temperature, coherent states remain coherent but their parameters 
evolve according to a{t) = a(0) e~3--*^*. A coherent superposition of two coherent states decohere as follows: 



\a{t)){f3{t)\ X e*''(*) e-'^^*) 



e{t) = (1 

d{t) = ^^(l-e--*) 



) 3(a/3) 



(71) 
(72) 

(73) 



The decoherence time is therefore given by a very simple formula: 

|a-/?|2 



'Dec 



7- 



(74) 



Using formula ^ in the Caldeira-Leggett case provides the same result and the temperature dependence is given, in 
the limit 7/2f7i ^ (f^i = ^f72 _^2/4) by. 



DT{t) = 



e(r)(i-e 



l + (e(T)-l).(l-e-^*) 



X 1 with e(r) = coth 



(75) 



Therefore, the t — > -\-oo limit of decoherence is independent of temperature but the decoherence time will scale with 
temperature according to the 6(T) factor: 



TDec(O) 



tanh 



2kBT 



(76) 



H. Decoherence at long time 

Applying previous results to the electromagnetic decoherence {I — 1 modes), we obtain: 

dl""Hn.n) = 4A,,,„.8in2(^:^) (77) 

As expected, the R/R limiting decoherence vanishes for ai = (T2 whereas the R/L one does not. The typical 
asymptotic value is proportional to A„^m. In fact, this number can be viewed as measuring the "distance" between 
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the two quantum states which built our Schrodinger cat. Using vertex operator with small values ofm and n in 
Schrodinger cats ^ and ^ produces mesoscopically separated coherent states in each mode. 

Figure 6 summarizes the electromagnetic decoherence exponent as a function of time for 7 1 ~ — 5 and for various 
temperatures. 

In the acoustic case (two dimensional phonon bath), contributions of all relevant modes should be summed. As 
before, the ai2 dependence disappears as soon as ai2 >> avs/cs- Introducing T = 2tt'^tlvs / {csgac{a)), the sum over 
all Luttinger modes up to the cut-off frequency can be evaluated: 

dWR)U) dWL)u) e-*/^-l , , 

1 + , (79) 



dWR) {00) dW^) {00) (t/T) 

where d'--^/-^^ {00) — d^^^^^^oo) = An^mCsL/{vsa). For the continuum approximation to be valid, we have assumed 
that L/a » vs/cs and therefore d(^/-^^(oo) >> 1. This also implies that most of the decoherence process is 
accomplished within the previously computed acoustic decoherence time 2T / d^-^^ {00) « T. 



V. CONCLUSION AND DISCUSSION 



Within the bosonization framework, we have shown how the coupling to an external quantum electromagnetic field 
or to a two or three dimensional bath of longitudinal phonons can lead to decoherence of a Schrodinger cat state 
formed of localized elementary excitations in a Luttinger liquid. Three different phases build the decoherence scenario 
for Schrodinger cat states of the Luttinger liquid (see 7 for the electromagnetic case) : 

• At a very short time, because of the time/energy uncertainty relation, energy exchanges between the environment 
and the Luttinger system are not conservative. In this regime, high frequencies take part in the decoherence 
process and non-Mar kovian effects are important. The precise time evolution of the system is strongly cut-off 
dependent. 

• After a transitory regime, the decoherence exponent becomes linear in time. This is the "Golden Rule" regime: 
energy conservation between the environment field and the Luttinger system is satisfied with a spectral width 
going down to its natural value (spontaneous photon or phonon emission). Only frequencies in resonance with 
the Luttinger eigenfrequencies contribute to dissipation. 

The structure of spectral densities for the environmental modes leads to a hierarchy of decoherence times 
corresponding to the multipolar expansion of the "radiation" emitted by the Luttinger system. In the acoustic 
case, the various decoherence times do not increase as much with increasing I. In the electromagnetic case, they 
increase with I. The total decoherence time is therefore much smaller in the acoustic than in the electromagnetic 
case. 

• At longer times, decoherence tends to saturation. In this regime, an effective Caldeira-Leggett model can be 
used to describe the dominating decoherence processes: in the electromagnetic case, one can keep only I = 1 
modes, corresponding to dipolar electromagnetic radiation. In the 2D acoustic case, one can use an effective 
Caldeira-Leggett model for all modes. Caldeira-Leggett computations are non-perturbative since they take into 
account all orders of the coupling between the quantum environment and the oscillator. 

The infinite time decoherence depends on two factors. The first one is, as expected, the distance between the 
two quantum states which depends on the dimension of the operators used in these states. The second one 
reflects the relative weight of each hydrodynamic mode in the decoherence process. In the electromagnetic 
case, since / = 1 modes dominate, it gives a geometrical factor depending on the relative position ai2. In the 
acoustic case, spatial dependence is almost always lost and we are left with an important mode number factor. 
That's why, although V'l?(o'i)|0) ^nd '0l?(ci)|O) can be considered as "mesoscopically close" with respect to their 
"distance" , decoherence is muclt^aster than dissipation in the acoustic case. This is a major difference with 
single mode decoherence studiesE3 where the ratio between decoherence and dissipation times is only due the 
distance between states entering the Schrodinger cat. 

Non-linearities in the spectrum of low energy excitations may also contribute to decoherence. As showed by 
HaldaneE3, non-linearities in the spectrum couple the bosonic modes of the theory, turning the simple free model 
used in bosonization into an interacting theory. Coupling between modes should also play an important role in 
the decoherence properties of Schrodinger cat states. Indeed, the model presented here provides an upper limit for 
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decoherence times. Non-linearities should be-also taken into account when investigating the resistive behavior of small 
metallic loops induced by inelastic collision^. 

In the present work, photons or longitudinal phonons initially at equilibrium were used as the thermal bath for 
the system, but of course, one could imagine various extensions. One could use another description for quantum 
fluctuations of the environment. For example, one may think about changing the state of the environment, taking for 
example into account an external microwave radiation. Increasing the incoming radiation power within the range of 
resonant frequencies should increase decoherence of Schrodinger cat states (enhancement of dissipation by stimulated 
emission of radiation). With such environmental states, one expects to meet also the problem of "decoherence 
repartition" between all the modes of the Luttinger system (even in the electromagnetic case) . Although this makes 
computations much harder to control, it may lead to more interesting behaviors. Finally, we are also investigating 
the extension of these ideas to two or three dimensional systems. 
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APPENDIX A: KELDYSH'S GREEN'S FUNCTIONS 

Let us recall the well known expressions for a single harmonic mode of frequency lu at temperature (3 (unit mass)@: 
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(Al) 

(A2) 

(A3) 

(A4) 
(A5) 

The retarded and advanced Green functions can be related to these expressions by 

Gn^Gr- G< and Ga ^ Gt ~ G> (A6) 
For the electromagnetic field, a mode decomposition can be used. Using results of appendix |^, we obtain: 

D^^^iix, t)- (y, s)) = Gu,At. s) (A7) 

The noise and dissipation kernels ^ and || can be expressed as: 

.-^..y) = coth(^) -°-(-;-(^^"--» (A8) 
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APPENDIX B: ELECTROMAGNETIC MODES IN A CYLINDRICAL CAVITY 

We shall decompose the transverse vector potential in orthonormal modes: 

A(r,t) = ^V^Kt)^/(r-) (Bl) 

These modes can be real, or complex. In the latter case, we assume the existence of an involution / i— > / over the set 
of indices implementing complex conjugation: (A/(r))* = Aj(r). Modes are normalized with respect of the volume 
V of the cavity by imposing an orthogonality condition: 

J Sr Ai{r).Aj{r)^5j j. 

The il}i{t) are coordinates for a set of independent harmonic oscillators of frequency w/. 

Expressions are given in cylindrical coordinates (r, {),z). R denotes the radius of the cavity. It is useful to introduce: 
A{r) ^ A{r)/s/V. 

TE modes Let I be the angular momentum of the mode around the Oz axis and kk its momentum along this axis. Let 
u'l „ denote the n-th zero of the derivative of the l-ih. order Bessel function Jj . Then, the orthoradial component 
of the TE [l, kz, n) complex mode is given by: 



We have lu 



TM modes Here u;^„ denotes the n-th zero of the l-ih order Bessel function J;. Then, the orthoradial component is 
given by: 



kzR 



(B3) 



We have uj — cyJ^^Tivd^^ZfR^ ■ 
Normalized real modes are inferred from these complex modes by: 



(+) 



V2 



{Ai + Aj) and A) 



(-) 



V2 



Ai) 



Expressions B2 and B3 show that: 



=Vicos{ld~kzz) 
A\~1 = -Vi sin {Id ~ kzz) 



(B4) 



(B5) 



where mathcalDj contains the r dependence and depends on the mode characteristics TE/TM and {l,kz,n). More 
precisely we have 



TE modes : Vi 
TM modes : P/ 



RkzV2 I JljjlUUn) 



(kzR) 



2 J'liuLn) 



(B6) 
(B7) 
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APPENDIX C: LONGITUDINAL ACOUSTIC MODES IN A CYLINDRICAL CAVITY 

The displacement field u{x, t) contains a gradient part: 

n(a=,i)=^^(V^.)(x) 

where the functions are eigenvalues of the Laplacian with von Neuman's boundary conditions. The modes ui = 
^ifi are normalized so that 



J (fr ui{r). uj{r) = 5j jV. 



The Qi{t) are coordinates for a set of independent harmonic oscillators with given frequency oj/. Longitudinal acoustic 
modes are indexed by / = {I, kz, n) for d = 3 and / = (/, n) for d = 2 and are typically of the form e*^'^~'^^^)A// where 

Mi contains the r dependence (w/ = cg + (uj nl^)"^ for d = 3 and = cgW; ^/R for d = 2): 

d = 2: Afr = -^^$4 (CI) 

d = S: Af, = "f"^ ^ff:'-"^ (C2) 



APPENDIX D: ASYMPTOTIC BEHAVIOR OF SPECTRAL DENSITIES 

Asymptotic expressions of the spectral densities corresponding to the infinite-cavity limit are useful to study the long 
time behavior of the decoherence and dissipation processes. The radial quantization is of Au!± = Au|^c/i? ~ ttc/R 
and the longitudinal one Au)\\ = 2Trc/h. We shall therefore use the following expressions for the extremas of the Bessel 
functions: 



Mui,n)-\ — — and J/(u;,„) ~ J (Dl) 



and the following low- and high- frequency expansions 

2H\ 



Mz):^^, and JK^)^ y^cos(^-/7r/2-7r/4) (D2) 



A straightforward algebra gives us: 

1. In the low-frequency regime {lotem << 1) 



JTE^^^ ^QED /U^Y (D4) 

'in c a \ 2iT J 

jnu;) = (D6) 
Henceforth, the total spectral densities behave as: 

" 3^ T ^ \-2^) 
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2. In the high-frequency regime [lotem >> 1) 



jn-) - ^ (D9) 



J-^^^H-w-i (DIO) 

These analytic results agree with the numerics, which are represented in figure 3. The main result is that all the 
modes are supraohmic at low frequency, expect the mode I = 1 which shows ohmic behavior. 

The same kind of expansion can be performed for the acoustic spectral densities. In this case, one obtains (d G 
{2,3}): 

• In the low frequency regime: 



Ji{uj)^gPh{L){uTsf'+'' (Dll) 



• In the high frequency regime: 



Ji{uj) ^ gpu{L) ^J-^ [uTsY-' (D12) 

Only the latter will be used to compute decoherence properties since uji falls into the high frequency regime since 
vs » cs- 

APPENDIX E: DECOHERENCE MATRIX COMPUTATIONS FOR A CALDEIRA-LEGGETT SOLUTION 

The appendix presents details of the computation of the decoherence matrix D ( ^ ) in the case of a damped oscillator 
solution. This has direct relevance for the Caldeira-Leggett model but also within ±hc_framework of the Breit-Wigner 
approximation in more general environments. These results have been discussed incj'EiJbut their derivatio n is re called 
here in a simpler way. Strictly speaking the approximate master equation approach described in section [VF is not 
valid for temperatures below the hj/ks- Therefore, this appendix aims at finding ultra- low temperature corrections 
to decoherence arising from the non-Markovian effects arising from the T — > behavior of the reservoir symmetric 
two point correlation function. 

Let fill and 7 denote the renormalized frequency of the oscillator and 7 the damping coefficient. We shall work in 
the weakly damped case, defined by — Qj^ — 7^/4 > 0, and measure the strength of dissipation by (j) such that 
tan (0) = 7/2r2jj. The decoherence matrix can then be expressed in terms of the following auxiliary functions: 

Z±{t,uj) = f e-7<^/2+'-(c^±f^i)'=rf<,^ (El) 
Jq 

S{t) - r(\Z+f + \Z^f){t,u;)J{uj)coth(^]duj, (E2) 



P{t) = / (Z+Z_)(i,c^)j7Hcoth('^ (E3) 







We have: 



1 f D12 



, = 4 fl + ,^^^f^ (s(^ + ^)+2'.X fpe-'.' (^,'«' -^)))+ (E5) 



-4PP) V \n nnj \ \n 

16cos2 ((Ale-^-y'X , , 
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and the coefRcients are given by: 

= " + c2^^7lo (^ + 2SR(Pe~^-("^*+^))), (E7) 

D12-D2, . _JfI!l__^(5smW + 23(Pe-*"''+*>)), (E8) 

The S and P functions can be computed by the residue theorem: 

S ^ - — ^- -T^ ^- + ; TT-, Jiiu) COth duj, ElO) 

1 /"^ 1 + e-''*e2*"i* - 2e-T*/2e*"i*e*'^* , , , //3fiw\ , , , 

P = 7^ -, 7^ ^^T7 7^ W)c0th -— duJ. (Ell) 

2 7_^ (cj + Oi+z7/2)(cj-^li-i7/2)^^ ^ V 2 / ^ ' 
The main contribution to S and P is due to the poles ±i7i + 17/2: 



= ^(l-e-^*)3fi(^:7(r!i+z7/2)coth(^^(r!i+i7/2) 



(E12) 
(E13) 

In this expression cut-ofF dependent quantities have been discarded since they can be shown to be of order 
(7/2171) log (A/J7/j), i.e. much smaller than the oscillator eigenfrequency's renormalization. The temperature de- 
pendent part also contain poles which give a series of exponentially decreasing terms of the form exp (— 27rni/rTft) 
(n > 1). At vanishing temperature, this series can be resummed into an algebraically decreasing correction as follows 
(here J{uj) — ^uj/it^I\): 

- -^^^5ft(e^"^*(zt5(zO+^"t5(-z-0) + e-^^^^*(z5(-zO+i"t5(i-,))), (E14) 

P^^°"-^ = - [S{zt)+S{-Zt) - S{zt) - S{-zt)) (E15) 

where 5(z) = Ei(l, 2) and z — (7/2 — ifJi) t. S{zt) has an asymptotic expansion in 1/t for t — > +00. 

At zero temperature, these corrective terms E14 and E15 dominate E12 and E13 for >> log(2ili/7). The very 
long time asymptotic of 5* is therefore given by: 

^c.Afi+40e-^y^cos(l)it)X ^^^^^ 



Qi V Trinity 

In case of weak damping 7 << fi/j, the following approximations can be made: (f> << 1, f2 ~ ~ fli and 
/(Oi + 17/2) ~ /(Sli). We will neglect P since, in full generality, P ^ 4)S. Plugging everything in D's expression, 
one finally ends up with a very long time asymptotics: 

D(T^Q)^[l-e +-e / j . ^ ^^^^ J . (E17) 

Henceforth, in the very weak damping situation, master equation results for the decoherence matrix can safely be 
extrapolated down to T = 0, even if strictly speaking non-Markovian effects should be taken into account. 
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Figure 1: The Keldysh contour consisting in a forward time path (upper branch) 
and a backward time bath (lower branch) 
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Figure 2: 

An annular ID Luttinger liquid embedded in a cylindrical electromagnetic cav- 
ity. 
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Figure 3: Normalized electromagnetic spectral densities for the 1 = 1,1 = 2 and 
I = 3 Luttinger modes (no cut off introduced here) . 
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Figure 4: Time dependence of the decoherence exponent d{t) and the renormal- 
ization coefRcient ri{t) for zero modes. 
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Figure 5: Spatial dependence of the electromagnetic decoherence time. The 
two top curves correspond to a T = 1 K whereas bottom ones correspond to 
T = 10 K. Here 6 = 27r(Ti2/L and fesTt = 2'KVsh/L. 
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